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Abstract 

Time evolution of the ionization state of metals in the cosmic baryons is investi- 
gated in a cosmological context without the assumption of ionization equilibrium. We 
find that a significant fraction of ionized oxygen (O Vll and O vill) in the warm-hot 
intergalactic medium (WHIM) is not in the ionization equilibrium state at a redshift 
of 2 ~ 0. We also investigate the effect on the detectability and observables of WHIM 
as a consequence of such deviation from ionization equilibrium. It is found that the 
detectability of WHIM is not altered very much both through its emission and ab- 
sorption signatures, but line ratios between O Vll and O Vlll are significantly different 
from those in the ionization equilibrium state. 
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1. INTRODUCTION 

After Fukugita, Hogan, & Peebles (1998) first pointed out the importance of studying 
the cosmic baryon budget, it is now widely accepted that only less than half of the cosmic 
baryons in the present universe have been detected, and the remaining portion, "the missing 
baryon" , has evaded the direct detection so far. Numerical simulations by Cen & Ostriker (1999) 
suggest that approximately 30 to 50 % of the cosmic baryons at 2 = are in the form of the 
diffuse intergalactic medium with temperature of 10^ K < T < 10^ K which is called warm-hot 
intergalactic medium (WHIM), and subsequent numerical simulations (e.g., Dave et al. 2001) 
with different numerical schemes and resolutions also support this picture consistently. Based on 
these numerical predictions, many observational efforts have been made to reveal the existence 
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of WHIM. 

Since the observational signature of WHIM is very weak, the detection of WHIM is 
observationally quite challenging and first proposed through its metal absorption features 
in the bright background beacons such as QSOs and blazers (Hellsten, Gnedin & Miralda- 
Escude 1998; Perna & Loeb 1998; Fang & Canizares 2000; Cen et al. 2001; Fang & Bryan 2001). 
After the first detection of Ovi absorption lines in the spectra of a bright QSO H1821+643 
by Tripp, Savage, & Jenkins (2000) and Tripp et al. (2001), number of detections are reported 
through absorption features of Ovi, Ovil, Ovill and Neix ions (Nicastro et al. 2002; Fang 
et al. 2002; Mathur et al. 2003; Fujimoto et al. 2004), but they are rather tentative. Only 
recently, the detection with sufficient signal-to-noise ratio is reported by Nicastro et al. (2005a) 
and Nicastro et al. (2005b) in which absorption signatures of WHIM are found at two redshift 
in the spectra of the blazer Mrk421 during its two outburst phases. Future proposed mis- 
sions such as Constellation- X and XEUS are expected to detect numerous WHIM absorbers. 
Detection of WHIM absorption in the spectra of afterglows of gamma-ray bursts (GRBs) are 
also proposed by Elvis & Fiore (2003) and Elvis, Nicastro & Fiore (2004) using a dedicated 
mission such as Pharos. Kawahara et al. (2005) investigate the feasibility of such detections in a 
realistic manner based on cosmological hydrodynamic simulations. As an alternative strategy, 
several tentative detections of WHIM through its metal line emission are claimed by Kaastra et 
al. (2003), and Finoguenov, Briel, & Henry (2003) using XMM- A/^ewton satellite. However, these 
detections are not significant enough to exclude the possibility that the observed emission lines 
are the Galactic ones because of the limited energy resolution (~ 80eV) of the current X-ray 
detectors. Recently, Yoshikawa et al. (2003), Yoshikawa et al. (2004), and Fang et al. (2005) 
show that future X-ray missions equipped with a high energy resolution spectrograph such as 
DIOS (Diffuse Intergalactic Oxygen Surveyor) and MBE (Missing Baryon Explorer) can de- 
tect the line emission of WHIM in a convincing manner. Future X-ray observations of WHIM 
through its emission and absorption features are expected to unveil the physical properties of 
WHIM beyond its existence. 

So far, most of theoretical and observational studies on WHIM assume that the cosmic 
baryons are in their ionization equilibrium under a given physical condition, irrespective of 
their thermal history. However, as shown in this paper, some portion of the cosmic baryons 
is expected to deviate from ionization equilibrium depending on its thermal history, since the 
timescales of ionization and recombination processes can be comparable or even longer than 
the age of the universe in a low density plasma such as WHIM. In this paper, we focus on the 
non-equilibrium ionization state of metals in WHIM. Since ionized metals such as O vi-vill 
and Neix are supposed to be good tracers of WHIM, the ionization state of such metals and 
their deviation from the ionization equilibrium state are of great importance in studying the 
detectability of WHIM and probing its physical condition. By relaxing the assumption that the 
cosmic baryons are in the ionization equilibrium, we follow the time evolution of the ionization 
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state of metals as well as hydrogen and helium in cosmological hydrodynamic simulations, 
and investigate to what extent the effect of the non-equilibrium ionization state affects the 
detectability and observational features of WHIM. 

The rest of the paper is organized as follows. The importance of non-equilibrium treat- 
ment of ionization balance in WHIM is described in section 2. In addition, technical details in 
computing time evolution of ionization fractions of ions are also shown in section 2. In section 
3, we consider a simple model of thermal histories of the cosmic baryons, and investigate the 
non-equilibrium ionization state of oxygen ions and their deviations from the ionization equi- 
librium state to understand general properties of the ionization history of oxygen in the cosmic 
baryons. The non-equilibrium ionization balance in the cosmological hydrodynamic simulation 
is shown in section 4, and the effect of the non-equilibrium ionization state on the detectabil- 
ity and observational features of WHIM is investigated in section 5. Finally, in section 6, we 
summarize our conclusion and discuss other important physical processes which may also affect 
the ionization state of WHIM. Throughout this paper, unless otherwise stated, we adopt a spa- 
tially flat cosmology with fim = 0.3, = 0.7, f2b = 0.04, and h = 0.7, where Q^n is the density 
parameter, Q\ the dimensionless cosmological constant, Qt, the baryon density parameter, and 
h the bubble constant in units of 100km s~^ Mpc~^. 
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Fig. 1. Redshift dependence of the cosmic UV and X-ray background radiation at photon energy of 0.5 
keV (left panel) and its spectrum at redshift z = 0, 0.5, 1, and 3 (right panel). In the right panel, we also 
show the observed background spectra in UV and X-ray bands obtained by ShuU et al. (1999) and Miyaji 
et al. (1998), respectively. 
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2. IONIZATION BALANCE 



The time evolution of ionization balance in coUisionally- and photo-ionized plasmas can 
be described by 

df- ^'"^ 

-i7 = Y.^i-k,kfk- Si-jjfj-ajfj + aj+ifj+i, (1) 

k=l i=j+l 

where j is the index of a particular ionization stage considered, Z the atomic number, fj the 
ionization fraction of an ion j. Sij denotes the ionization rate for an ion j with the ejection 
of i electrons, and aj is the recombination rate of ion j. As ionization processes, collisional. 
Auger, charge-transfer and photo-ionizations are taken into account. Recombination processes 
are composed of radiative and dielectric recombinations. The ionization equilibrium state, j'P, 
can be obtained by solving 

j-i z+i 

k=l «=i+i 

Ionization and recombination rates adopted in this paper are calculated by utilizing the SPEX^ 
ver 1.10 software package. 

As a source of photoionization, we adopt UV and X-ray background radiation calculated 
using the CUBA code^ (Haardt & Madau 2001) throughout this paper. The spectra and redshift 
evolution of the adopted background radiation are shown in figure 1. It should be noted that the 
intensity of the background radiation rapidly drops at redshift z <2, while the overall spectral 
shapes do not change so much. In the right panel, we overlay the observed UV and X-ray 
background radiation estimated by ShuU et al. (1999) and Miyaji et al. (1998), respectively, 
for comparison. It can be seen that the spectrum of adopted background radiation at 2; = 
is consistent with the observed ones within a factor of two. Although ionization fractions of 
ions in the collisional ionization equilibrium entirely depend on gas temperature, they depend 
on gas density as well as gas temperature under the presence of photoionizing background 
radiation, and the effect of photoionization by the UV and X-ray background is important in a 
low density and low temperature plasma. Figure 2 shows the ionization fractions of O vi, O vil 
and Ovill ions in collisional ionization equilibrium (i.e. no background radiation) and those 
in ionization equilibrium under the presence of the UV and X-ray background at 2; = 0. It is 
obvious that the effect of photoionization is quite important at a typical temperature range of 
WHIM (10^ K < T < 10^ K) and that it increases the ionization fractions of ionized oxygen in a 
lower density regime. 

According to equation (1), the timescales of ionization and recombination for an ion j, 
Tion and Tree respectively, can be estimated by 

^ http://'W'W'w. sron.nl/divisions/hea/spex/ 

^ http://pitto.mib.infn.it/haardt/cosmology.html 
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Fig. 2. Ionization fractions of O vi, O vii and O viii in ionization equilibrium as a function of temperature 
for the low density plasma with tih ~ 10^^ and 10^^ cm^"^ under the presence of photoionizing background 
radiation at a redshift of z = 0. Ionization fractions in coUisional ionization equilibrium are also shown for 
the comparison. 




(3) 



and 



Tree = a/- (4) 

Figure 3 shows temperature dependence of Tjon and Tree of O vi, O Vll and O Vlll ions for 
hydrogen number density of nn = 10"^ and 10"^ cm~'^ under the presence of the photoionizing 
background radiation at a redshift of 2; = 0. In a temperature range T ^ 10^ K, the recombination 
timescales are longer than the ionization timescales, and are comparable to or even longer than 
the bubble time, H^^ (~ 13 Gyr), where = lOO/i kms~^ Mpc~^ is the bubble constant. In a 
lower temperature regime with T ^ 10® K, the recombination timescales are still comparable to 
the bubble time for hydrogen number density of rin = 10~® cm~^. The ionization timescales in 
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the lowest temperature range (~ 10^ K) are determined by the intensity of the UV and X-ray 
background radiation. It can be seen that the ionization timescales of O vil and O vill ions are 
also comparable to the bubble time at a temperature of T ^ 10^ K. Actually, these behaviors 
of ionization and recombination timescales are also the case for hydrogen- and helium-like ions 
of other metal such as carbon and nitrogen. Considering these facts, one can easily expect 
that the ionization balance is not in the equilibrium state in a low density plasma like WHIM. 
Therefore, for proper treatment of ionization balance in WHIM, we have to follow the time 
evolution of ionization fractions by directly integrating equation (1). 
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Fig. 3. Timescales of ionization and recombination for Ovi, Ovn and Ovni ions in hydrogen number 
density of nn = lO^^cm^"^ (left) and lO^^cm^"^ (right) under the presence of UV and X-ray background 
radiation at a redshift of z = 0. Horizontal dotted lines indicate the hubble time H^^ (~ 13Gyr). 

Since the timescales for the ionization and recombination differ by many order of mag- 
nitude depending on ionic species and ionization stages, equation (1) is a stiff set of differential 
equations. In numerically solving a stiff set of equations, explicit integration schemes are 
unstable unless unreasonably small timestep is taken, which leads to prohibitively expensive 
computations. Therefore, implicit schemes are required to numerically integrate equation (1). 
In this paper, we adopt a backward difference formula (BDF), which is also adopted in Anninos 
et al. (1997) for solving chemical reactions of atomic and molecular hydrogen in the primordial 
universe. In the BDF, ionization fractions of all ionic species are sequentially updated in the 
order of increasing ionization states rather than being updated simultaneously, and the source 
terms contributed by the ionization from and recombination to the lower states are evaluated 
at the advanced timestep. The discretization of equation (1) can be written as 



t+At 



J 



(5) 



We confirm that this scheme is stable and as accurate as the eigenvalue decomposition scheme by 
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Kaastra & Jansen (1993) adopted for solving the non-equilibrium ionization state in supernova 
remnants, but our scheme is found to be computationally less expensive. In this paper, we 
simultaneously solve the time evolution of ionization fractions of H, He, C, N, O, Ne, Mg, Si, 
and Fe, but we focus on only the ionization states of oxygen. 

Throughout this paper, we assume spatially uniform and time independent metallicity 
of Z = 0.1^0, where Zq denotes the solar abundance. Note that the metallicity Z is very 
insensitive to the resulting non-equilibrium ionization state as long as Z <C 1. Of course, in 
reality, metals in intergalactic medium are ejected from galaxies through the feedback of star 
formation activities, and thus the distribution of metals is inhomogeneous and time dependent. 
Actually, for a fully consistent treatment, such feedback effect must be included in the rhs 
of equation (1) as a source term. In this paper, however, since we are interested in the non- 
equilibrium ionization state that arise along the thermal histories of baryons in a context of 
cosmological structure formation, we adopt such an assumption and believe that it is acceptable 
for qualitative studies of ionization states of WHIM. 

3. SIMPLE MODELS 

In this section, we explore the general behaviors of the non-equilibrium ionization state 
using simple thermal histories of the cosmic baryons. Here, we consider a series of thermal his- 
tories depicted in figure 4, in which baryonic matter has hydrogen number density nH,o(l + ^1)^ 
and temperature T, at an initial redshift of Zi = 10, keeps its comoving density and temperature 
constant before increasing its density and temperature to riu = nH,o(l + ^s) and T = Tg, respec- 
tively, by experiencing shock heating at redshift z = Zg, and keeps the postshock thermal state 
till z = 0. A physical scenario we have in mind is that the baryonic matter is diluted according 
to the cosmological expansion until it experiences the cosmological shock heating, and that it 
gets quasi-virialized inside collapsed halos or filamentary structures after the shock heating. 
Here ?7,h,o is the present-day hydrogen number density, and we set the initial temperature to 
Tj = 5 X 10'^ K. Therefore, these thermal histories can be parametrized by three parameters, the 
postshock overdensity 6s, postshock temperature Tg, and the redshift of shock heating Zg. In a 
cosmological model adopted in this paper, we have ?t,h,o = 2.2 x 10"^ cm~'^. We assume the UV 
and X-ray background radiation described in the previous section as a photoionization source, 
irrespective of the adopted thermal histories. 

Actually, the thermal history of the cosmic baryon in a cosmological hydrodynamic 
simulation is fairly well represented by this simple model. Figure 5 shows the time evolution of 
temperature and hydrogen number density of gas particles from a cosmological hydrodynamic 
simulation presented in section 4. In both panels, time evolutions of five randomly selected 
gas particles inside a filamentary structure at z = are shown. One can see that each gas 
particle have nearly constant temperature and physical density proportional to (1 + z)~^ before 
it experiences the first shock heating. After the first shock heating, although thermal histories 
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of gas paricles in the simulation is not so simple as that presented in figure 4, it can be a 
viable model for qualitative description of the thermal history of the cosmic baryon since the 
timescales for temporal variations in density and temperature are much shorter than Tjon and 
Tree of oxygeu ious (O VI, O VII and O Vlll). It should be noted that, in figure 5, we can see that 
some particles experience multiple shocks in a sufficiently short timescale compared with Tjon 
and Tree, and increase their postshock density by a factor larger than 4, the value in the strong 
shock limit. Therefore, we allow the density jump a.t z = Zs to be greater than a factor of 4 in 
order to bracket the plausible range of the thermal histories of the cosmic baryons. 
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Fig. 4. Time evolution of density and temperature of baryons for the model calculations. We have three 
free parameters Ss, Tg and Zg. (see text for detail.) 

At the initial redshift Zi = 10, we assume that all baryons are in the ionization equilibrium 
state. This assumption can be justified by the fact that physical gas density is so high at such 
a high redshift that ionization equilibrium can be quickly realized. Equation (1) is integrated 
along the thermal histories described above till z = for a given set of parameters 6s, Tg, and 
^s- Figure 6 shows the ratios of the resulting ionization fractions at 2; = for thermal histories 
with Zs = 1.0 relative to those in the ionization equilibrium state, fj/f'P, for O vi, O vil, O vill 
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Fig. 5. Time evolution of temperature (left panel) and hydrogen number density (right panel) of gas 
particles in a cosmological hydrodynamic simulation presented in section 4. See section 4 for details of the 
simulation. 

and O IX on a (1 + 6s)-Ts plane. White contour lines are separated by a factor of 1.1 in ratios 
of ionization fractions, fj/fj'^, where the bold lines are contours for fj/f'P = 1- Black contours 
with numeric labels indicate the ionization fractions in the ionization equilibrium state on 
{1 + 6s) — T plane. One can see that, for O vi ion, the ratios of ionization fractions significantly 
deviate from unity only in the region where Ovi fraction in the equilibrium state, /ovi; is 
negligible, and close to the equilibrium value in the region where /qvi is higher than 0.1. This 
is because, as can be seen in figure 3, timescales for both of ionization and recombination of 
O VI ions are much shorter than the bubble time in the density and temperature ranges where 
/ovi is significant. On the other hand, for Ovil, Ovill and Oix, we have significant deviation 
from the ionization equilibrium state in those density and temperature ranges in which their 
ionization fractions are not negligible, especially at low density regions with 1 + (5s < 1.5. Only 
the region with 1 + 5s > 10 and Ts > 10^'^K is almost free from significant deviation from the 
ionization equilibrium state. 

Figure 7 shows the time evolution of ionization fractions of O vil, O vill and O ix for 
thermal histories with four sets of 1 + 5s and Tg marked in figure 6 (A, B, C and D) and Zs = 1.0. 
For all the thermal histories, ionization fractions of O Vlll increase rapidly just after the shock 
heating, while those in the ionization equilibrium state are lower than the preshock ionization 
fractions. This is because the timescale for coUisional ionization to O vill ion is much shorter 
than that of recombination timescale. After taking their maximum values, O vill fractions 
decrease toward the equilibrium fractions. Ionization to O ix ions slightly slower and is delayed 
compared with that to O Vlll ions. Ionization fractions in the equilibrium state also change 
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in time because the intensity of the UV and X-ray background radiation is decreasing at a 
redshift z <2. In the history A with 1 + 5s = 10 and T = lO*" K, O vill and Oix fractions 
in the ionization equihbrium state are increasing and decreasing in time, respectively. In the 
non-equilibrium ionization state, however, the ionization fractions of these ions cannot catch 
up with the ionization equilibrium state because of slow recombination processes from O ix to 
O VIII. Therefore, it can be said that the decrease of the UV and X-ray background radiation 
at 2; ^ 2 is so rapid that it eventually enlarge the deviation from the ionization equilibrium 
state. The similar behavior can be also seen in the histories C and D. On the contrary, in the 
history B with higher density and temperature than A, the ionization equilibrium is achieved 
because ionization processes are dominated by collisional ionization and the decrease in UV and 
X-ray background intensity does not significantly change the equilibrium ionization fractions. 
Figures 8 and 9 are the same as figure 6 but for results with = 0.5 and 0.1, respectively. 
As is the case for the results with Zg = 1.0, we also have significant deviation from ionization 
equilibrium for O Vll, O Vlll and O ix ions, while O vi fractions are close to the ionization 
equilibrium states for the region where /qvi > 0.1. For Zs = 0.5 and 0.1, fractions of O Vll and 

VIII are larger than the equilibrium fractions, while O ix ion is under-ionized in regions with 

1 + 5s < 10 and Tg > lO^K. These features can be understood by the same argument as the result 
for a history with 2s = 1.0 that the timescale for ionization to O vill is much shorter than that 
of recombination from O VIII to O Vll, and that ionization to O ix is slower than that to O Vlll. 

4. RESULTS IN COSMOLOGICAL SIMULATION 

In this section, we compute the non-equilibrium ionization state for more realistic ther- 
mal histories of baryons obtained from a cosmological hydrodynamic simulation. The simulation 
data of Yoshikawa et al. (2003) is used in this paper, which is computed using a hybrid code 
of a Particle-Particle-Particle-Mesh (PPPM) poisson solver and smoothed particle hydrody- 
namics (SPH) scheme. The details of our code are described in Yoshikawa et al. (2001). We 
adopt a spatially flat ACDM universe with = 0.3, = 0.7, fib = 0.015/i~^, as = 1.0 and 
h = 0.7, where as is the r.m.s. density fluctuation smoothed over a scale of 8/i~^Mpc. We 
employ 128'^ particles each for dark matter and gas particles within a periodic simulation box 
of 75h~^Mpc per side. Therefore, the mass of each dark matter and gas particle is 2.2 x lO^^M© 
and 2.4 x lO^M©, respectively. The effect of radiative cooling is incorporated using the cooling 
rate for a metallicity of [Fe/H]= —0.5 by Sutherland & Dopita (1993). The effect of energy 
feedback from supernovae is ignored in our simulation. 

The initial condition is setup at z = 36 and we have 51 outputs from z = 9 to z = with 
an equal interval in ln(l + z). For each gas particle, we integrate equation (1) till z = using the 
BDF described by equation (5) using its thermal history, or more specifically, a history of its 
density and temperature. As an intial condition of equation (1), the output at 2 = 9 is adopted. 
As in the simple model, it is assumed that all baryons are in the ionization equilibrium state 
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Fig. 6. Maps of ratios of ionization fractions between non-equilibrium and equilibrium ionization states 
for O VI, Ovil, Ovili and Oix at redshift z = on a (1 + 5s)-Ts plane, where shock heating epoch is set 
to Zs — 1.0. White contours are separated by a factor of 1.1 in the ratios of the ionization fractions. Black 
lines with numeric labels are the contours of the ionization fractions in the ionization equilibrium states. 

in the initial epoch. Thermal history of each gas particle is interpolated if an time interval 
between two continuous outputs is longer than the timestep At in integrating equation (5). 
The calculations are carried out with timesteps of 10^ yr and 3 x 10^ yr for a convergence 
test, and we verify that the results with two different timesteps are consistent with each other. 
Note that our methodology is not fully self-consistent because the thermal history of each gas 
particle, with which its non-equilibrium ionization state is calculated, is computed with the 
cooling rate in the coUisional ionization equilibrium by Sutherland & Dopita (1993). However, 
since the cooling time for hot (T ~ lO^K) and low density (1 + 5b — 1) gas like WHIM is much 
longer than the bubble time, the effect of radiative cooling on its thermal history is rather mild, 
and such inconsistency in radiative cooling does not drastically change its thermal history. 
Therefore, we believe that our results presented in this paper is useful as a first approximation 
of non-equilibrium ionization state in WHIM. 

Figure 10 depicts mean ratios of ionization fractions of O vi, O Vll, O Vlll and O ix in the 
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Fig. 7. Time evolution of ionization fractions of Ovii, Oviii and O ix for four sets of the densities and 
temperatures at redshift z = marked in figure 6. Bold and thin lines indicate the non-equilibrium and 
equilibrium states, respectively. 

non-equilibrium state at 2; = relative to those in the equilibrium state on a (1 + 6h)-T plane, 
where is the baryon overdensity. Black contours in each panel indicate the mass distribution 
of the corresponding ion, and 25%, 50% and 75% of its mass are enclosed by three contour 
lines from inside to outside. One can see that most of Ovi ions are close to their ionization 
equilibrium states, and that, the deviations of Ovil, Ovill and Oix ions from the ionization 
equilibrium are not negligible. In addition, all of O vil, O Vlll and Oix ions with (1 +5b) > 10^ 
and T > are close to the ionization equilibrium. These features are similar to the results 

for simple thermal histories described in the previous section, and can be understood by the 
same arguments on the timescales of ionization and recombination processes. It should be also 
noted that deviation from ionization equilibrium states at low density regimes with l + 5b< 10 
is quite similar to that in the simple models with = 0.1 (see figure 9), indicating that the 
baryons in these low density regions recently experienced shock heating. 

In figure 11, we show the spatial distribution of hydrogen number density and the ratios 
of O VI, O VII and O vill fractions in the non-equilibrium state relative to those in the ionization 
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Fig. 8. Same as figure 6 except that shock heating epoch is set to Zg = 0.5. 

equilibrium state in a 7.5h~^ thick slice taken from the simulation box at 2; = 0. In this slice, 
we have a galaxy cluster at {X,Y) = (32, 25)/i~^ Mpc, and a prominent filamentary structure 
associated with this galaxy cluster. A close look at figure 1 1 shows that we have significant de- 
viation from the ionization equilibrium state at the outskirts of galaxy clusters and filamentary 
structures, where cosmological shocks occur (Kang et al. 2005). 

5. EFFECT ON DETECTABILITY AND OBSERVABLES OF WHIM 

The deviations of ionization fractions of O vi, O vil and O Vlll from the ionization equi- 
librium state shown in the previous section will affect the observational features of WHIM. 
In this section, using the results obtained in the previous section, we present to what extent 
the detectability and observables of WHIM through its emission and absorption signatures are 
altered by relaxing the assumption of the ionization equilibrium. For a mock observation of 
WHIM, we construct a lightcone data up to a redshift of z = 0.3 by stacking 11 simulation 
cubes in the same way as described in Yoshikawa et al. (2003). The lightcone data has a 5° x 5° 
field-of-view and contains the ionization fractions of ions for each gas particle computed in the 
previous section as well as those in the ionization equilibrium state. 
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Fig. 9. Same as figure 6 except that shock heating epoch is set to Zg = 0.1. Note that the ranges of the 
ratios are different from tliose in figure 6 and 8. 

5.1. Detection through emission 

Surface brightness of emission lines of Ovil at 574eV and Ovill at 653eV (^ovii and 
'S'oviii! respectively) are computed for both of the equilibrium and non-equilibrium states. 
Here, we assume the latest specification of the detector proposed for DIOS mission (Ohashi 
et al. 2004), namely a spectroscopic energy resolution of 2 eV and a spatial resolution of 5 
arcmin. Figure 12 depicts the ratios of surface brightness in the non-equilibrium state relative 
to that in the ionization equilibrium state 5'ovii/5'ovn and 5'oviii/>S'ovjjj, separately for line 
emitters at redshift ranges of 0.0 < z < 0.15 (lower panels) and 0.15 < z < 0.3 (upper panels). 
Vertical dotted lines in all the panels indicate the nominal detection limit of DIOS mission 
(10~^^ ergs~^ cm~^ sr~^). Note that we have significant deviations of surface brightness ratios 
from unity at ~ 10~^ ergs~^ cm~^ sr~^ in O VII emission and ~ 10"*^ ergs~^ cm~^ sr~^ in O VIII 
emission. It is found that these strong deviations take place at the very central regions of 
galaxy clusters and are caused by shocks inside them. One can see that surface brightness of 
Ovill (653 eV) emission lines in the non-equilibrium ionization state is quite close to that in 
the equilibrium state, and that O vil (574 eV) emissions are only slightly brighter than those 
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Fig. 10. Mean ratios of ionization fractions at 2 = between non-equilibrium and equilibrium ionization 
states for O vi, O vii, O viii and Oix on the (1 + 5t,)^T plane. White contours are drawn in the regions 
with 0.5 < fj/fj'^ < 2 and separated by a factor of 1.2 in the ratios of ionization fractions. The white solid 
lines indicate the contour for fj/fj'^ = 1. Three black contours from inside to outside enclose the 25%, 
50% and 75% of mass of each ionic species. 

in the ionization equilibrium state above the detection hmit of BIOS mission, irrespective of 
redshift of the hne emitters. Therefore, the detectabihty of WHIM through its oxygen emission 
hues is not significantly different from the theoretical predictions based on the assumption 
of ionization equilibrium. Actually, the previous works on oxygen line emission of WHIM by 
Yoshikawa et al. (2003), Yoshikawa et al. (2004) and Fang et al. (2005) all assume the collisional 
ionization equilibrium. Therefore, it is useful to understand to what extent such an assumption 
is acceptable. Figure 13 shows the comparison of the surface brightness of Ovil and Ovill 
emissions in the non-equilibrium state with those in the collisional ionization equilibrium, and 
indicate that the assumption of collisional ionization equilibrium is fairly valid for oxygen line 
emitters brighter than the nominal detection limit of BIOS, 10~^^ erg s~^ cm~^ sr~^. Therefore, 
the numerical predictions of detectabihty of WHIM through its oxygen line emission presented 
in Yoshikawa et al. (2003) and Yoshikawa et al. (2004) is not significantly altered by the effect 
of the non-equilibrium ionization balance. 
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Fig. 11. Spatial distribution of hydrogen number density, and ratios of ionization fractions of O VI, O Vii 
and O VIII between non-equilibrium and equilibrium states in the simulation volume at z = 0. The thickness 
of the slice is 7.5/i~^Mpc. Contours of baryon mass density for 1 + i5b = 1 are also shown in the top-right, 
bottom- left and bottom right panels. 

5.2. Detection through absorption 

Column densities of Ovi, Ovil and Ovill {Nqyi, Nqvu, and A^'oviii, respectively) and 
those in ionization equilibrium (iVovi, ^ovii ^oviii) calculated for 10^ randomly se- 
lected line-of-sights in the 5° x 5° field-of-view by computing line integrals through the smooth- 
ing kernel of each gas particle, and individual absorbers are identified in the same way as done 
in Chen et al. (2003). Each line-of-sight is divided into ~ 16000 bins along redshift direction. 
Thus, each bin has a width of ~ 5 km s~^ in velocity units. Figure 14 shows the column 
density ratios of the identified absorbers in the non-equilibrium state relative to those in the 
ionization equilibrium state (iVovi/-^ovi' ^ovii/^ovii ^oviii/^^oviii) ^ function of the 
latter for absorbers at redshift ranges of 0.0 < z < 0.15 (lower panels) and 0.15 < 2 < 0.3 (upper 
panels) separately. As expected from the results shown in previous sections that most of O VI 
ions are quite close to the ionization equilibrium state, column densities of O VI absorbers are 
almost the same as those in the ionization equilibrium state, irrespective of the redshift of the 
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Fig. 12. Ratios of surface brightness of Ovii (574eV) and Oviii (653eV) emission lines for 
non-equilibrium states relative to equilibrium states as a function of surface brightness in equilibrium 
states. Upper and lower panels are for line emissions at redshift ranges of 0.15 < z < 0.3 and 0.0 < z < 0.15, 
respectively. Symbols and error bars indicate the averages and the standard deviations for each horizontal 
bin. Vertical dotted lines in both panels indicate the nominal detection limit of DIOS mission. 

absorbers. On the other hand, we have systematically larger O vil and O vill column densities 
in the non-equilibrium ionization state and also have large statistical dispersions. However, the 
extent of difference in column densities between equilibrium and non-equilibrium states is so 
mild that the effect on detectability of WHIM through oxygen absorption lines is negligible. 

5.3. Line Ratios 

Apart from intensities of emission and absorption lines of WHIM, their ratios are also 
a very important observables as diagnostics for the nature of WHIM, in particular its tem- 
perature. Deviation from ionization equilibrium will cause difference in ratios of emission and 
absorption line intensities of WHIM from those with the assumption of ionization equilibrium, 
and uncertainty in estimation of temperature of WHIM. 
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Fig. 13. Same as figure 12, except that ratios relative to surface briglitness in the coUisional ionization 
equilibrium state are shown. 

Figure 15 depicts the relation of emission line ratios S'oviii/'S'ovii between equilibrium 
and non-equilibrium states for line emitters at 0.0 < z < 0.15 (left panels) and 0.15 < z < 0.3 
(right panels). Here, line emitters with either S'qvh or 5'ovin larger than 10~^^ ergs"^ cm~^ sr~^ 
are considered. One can see that the line ratios 5'oviii/5'ovii in the non-equilibrium state are 
systematically lower than those in the ionization equilibrium state, especially at S'oviii/'S'ovii ^ 
10 for both redshift ranges. Actually, we find that line emitters with Soviii or Sovii — 10"^^ 
ergs~^ cm~^ sr~^ mainly contributes to these differences line ratios, and that line ratios of 
emitters whose Sovu and Soviii exceed 10~^° ergs~^ cm~^sr^^ are statistically consistent with 
those in the ionization equilibrium state. 

Relation between column density ratios Nqyu/ Nqyi and A^^oviiiZ-^ovii in the ionization 
equilibrium and non-equilibrium states are shown in figure 16 for absorbers at redshift 0.0 < 
2 < 0.15 (left panel) and 0.15 < 2; < 0.3 (right panel), separately. Here, column density ratios 
Nqyu/Novi (Novm/Novn) are considered for absorbers whose A^ovi or A^ovii (A'oviii or A^'ovii) 
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Fig. 14. Ratios of Ovi, Ovii and Oviii column densities of absorbers at redshift 0.0 < z < 0.15 (lower 
panels) and 0.15 < z < 0.3 (upper panels) in non-equilibrium states relative to those in ionization equilibrium 
states identified in the lightcone data as a function of column densities in ionization equilibrium states. 
Symbols and error bars indicate the averages and standard deviations for corresponding horizontal bins, 
respectively. 
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Fig. 16. Relations of column density ratios A^ovii/-^ovi (filled squares) and A^oviii/-^OVii (open trian- 
gles) between non-equilibrium and equilibrium states. Symbols and error bars indicate the averages and 
standard deviations for corresponding horizontal bins, respectively. 

exceed 10^'' cm~^. Although Nqvu/ Nqyi in the non-equihbrium state statistically agrees well 
with that in the ionization equilibrium state in both redshift ranges, the non-equilibrium state 
exhibits systematically lower A^oviii/^^ovii than the ionization equilibrium state. It is found 
that these differences in iVoviii/A''ovii are due to the absorbers with low (~ lO^^cm"^) column 
density of O vil and O vill ions. 

The location and the physical conditions where the line ratio between O Vlll and O Vll 
deviates from that in ionization equilibrium state are of great interest. Here we define the ratio 
of ionization fractions of O vill relative to O vil in equilibrium and non-equilibrium state as 

Jovii 



and 



cq 

i?^^ = (7) 

Note that 5'oviii/5'ovii and Noym/Noyu are roughly proportional to R. Figure 17 shows the 
map of R/W^ at a redshift of z = in the same slice of the simulation volume as figure 11, and 
indicates that we have R/R'^^ ^ 1 ubiquitously in filamentary structures. This is why the line 
ratios, 5'oviii/5'ovii and A^'oviii/A'^ovii! are systematically smaller than those in the ionization 
equilibrium state. The left panel of figure 18 shows the contour map of R/R'^'^ on a (1 + 5b)^ 
T plane, on which three contours for i?®^ = 0.1, 1 and 10 are overlayed from bottom to top. 
Considering the fact that the 5'oviii/5'ovii and A^'oviii/A'^ovii in the non-equilibrium ionization 
state are lower than their equilibrium values at Soyjij/Sovu ^ 10 and iVoviii/^ovii ^ 
spectively, the left panel of figure 18 indicates that the regions responsible for such departures 
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Fig. 17. Spatial distribution of R/ R^^ at z = in the same slice of the simulation volume as in figure 11. 
Contours for baryon mass density of 1 + (5b = 1 arc also shown. 




Fig. 18. Maps of R/R'^^ at z = on a (1 +6h)-T plane in a numerical simulation (left panel) and in a 
simple model with Zg = 0.1 (right panel). White contours are separated by a difference of 0.25 in R/R'^^ 
and bold lines indicate R/R'^^ — 1. Three black contours in each panel are for R'^'^ = 0.1, 1 and 10 from 
bottom to top. 
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in line ratios from the equilibrium state have a typical density of 1 + 5b — 10 and temperature 
of T ~ lO^^^K. The right panel of figure 18 is the map of R/ R^"^ at 2 = in the simple model 
with = 0.1 described in section 3, and shows the essentially same behavior of R/R^'^ on the 
(1 + 5\^)-T plane. This suggests that the regions with R/ R^^ < 1 experienced a shock heating 
very recently around a redshift of z ~ 0.1, and that over abundant O VII ions due to the non- 
negligible timescale for the ionization from O vil to O Vlll reduces R/R°^ to below unity. Note 
that the loci with R/ R"^"^ < 1 on the (1 + 6h)-T plane in figure 18 and those with /ovii/ /ovii > 1 
roughly coincide (see upper right panels of figure 9 and 10). 

6. CONCLUSION AND DISCUSSION 

We investigate the non-equilibrium ionization state of oxygen ions in WHIM by relaxing 
the assumption of ionization equilibrium and directly integrating the time evolution of their 
ionization fractions along the cosmological thermal histories of the baryons. 

First, we consider the general features of the non-equilibrium ionization state for simple 
models of thermal histories of the cosmic baryons, and find that ionization states of O vil, 
O VIII and O ix ions are significantly different from those in the ionization equilibrium state 
especially in typical physical conditions for WHIM, while most of O vi ion are close to the 
ionization equilibrium state. These deviations from the ionization equilibrium are caused by 
the fact that timescales for ionization and recombination processes of these ions are comparable 
to or even longer than the hubble time in low density plasma like WHIM. Furthermore, we find 
that rapid decrease of the intensity of the UV and X-ray background radiation at redshift z <2 
also enlarge the deviation from the ionization equilibrium state. This is because the timescale 
of the recombination processes is so long that the recombination processes cannot catch up 
with the ionization equilibrium state, in which ionization fractions also vary with time as the 
UV and X-ray background radiation decays. 

The evolution of ionization fractions of oxygen ions for thermal histories obtained in a 
cosmological hydrodynamic simulation is also investigated. The ionization state of oxygen ions 
at redshift of 2 = is similar to that obtained in the simple models of thermal histories in which 
cosmological shocks take place at relatively recent epoch, say z ^ 0.5. Statistically, O Vll and 
O VIII ions are over-abundant compared with the ionization equilibrium state in typical physical 
conditions for WHIM, while Ovi ions are nearly in the ionization equilibrium state. Spatial 
distribution of the oxygen ions reveals that the deviations of ionization fractions of O vil and 
O VIII ions from the ionization equilibrium state take place at the outskirts of galaxy clusters 
and the edge of filamentary structures. 

We also investigate the effect of such deviations of oxygen ions from the ionization equi- 
librium state on the detectability and observables of WHIM through its emission and absorption 
line features. It is found that although the surface brightness and column density of oxygen 
emission and absorption lines slightly differ from those in the ionization equilibrium state, the 
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overall detectability of WHIM is almost the same as what we have under the assumption of the 
ionization equilibrium. However, the ratios of emission and absorption line intensities between 
O VII and O vill, which are considered to be important probes of temperature of WHIM, are 
significantly different from those in the ionization equilibrium. This means that the observed 
line ratios of WHIM may not reflect its real temperature. It is found that the regions where the 
line ratios are different from those in the ionization equilibrium are located at outer edges of 
filamentary structures. Distribution of such regions on a density-temperature plane is similar 
to the result of the simple model for the cosmological thermal history with = 0.1 described 
in section 3. These facts indicate that recent cosmological shocks near the edge of filaments 
around a redshift of 2 ^ 0.1 cause the departure of line ratios from the ionization equilibrium 
state. 

In this paper, several important physical processes are missed. First of all, in our simula- 
tion, the effects of energy feedback and metal ejection from galaxies are ignored, and we assume 
that the metallicity of the cosmic baryons is spatially uniform and constant irrespective of red- 
shift. Since oxygen atoms are mainly supplied by type-II supernovae followed by the formation 
of massive stars, the spatial distribution of oxygen is obviously time dependent and inhomoge- 
neous. Furthermore, the energy feedback by supernova explosions will affect the thermal his- 
tory of intergalactic medium (or WHIM) in the vicinity of galaxies. Therefore, a self-consistent 
treatment of galaxy formation and its feedback effect is required for better understanding of 
the non-equilibrium ionization state of ions in WHIM. Recently, Cen & Ostriker (2006) and 
Cen & Fang (2006) performed numerical simulations in which the non-equilibrium ionization 
state of oxygen ions is considered together with energy feedback and metal pollution caused 
by galaxy formation. According to their results, energy feedback induces strong shocks and 
additional non-equilibrium ionization state in the intergalactic medium around galaxies if the 
energy feedback is strong enough as performed in their simulations. Therefore, the deviation 
from ionization equilibrium of oxygen ions in WHIM presented in this paper can be regarded as 
its lower limit, and the actual deviations from the ionization equilibrium may be much larger. 

Secondly, it is assumed that equipartition between electrons and ions is achieved instan- 
taneously, or, in other words, electrons and ions always have a locally common temperature. 
At cosmological shocks, by which most of WHIM is heated, nearly all of the kinetic energy is 
converted into the thermal energy of ions just after the shock. Actually, collisional relaxation 
between ions and electrons through coulomb scattering takes a long time which is comparable 
to the bubble time in a low density plasma like WHIM. Therefore, as numerically shown by 
Takizawa (1998) and Yoshida, Furlanetto, & Hernquist (2005), the electron temperature Tg 
can be systematically lower than the ion temperature T; by a factor of a few in the outskirts 
of galaxy clusters and WHIM, unless other physical processes which convert the bulk kinetic 
energy into electron thermal energy (e.g.. Laming 2000) work efficiently. Since the ionization 
and recombination rates Sij and aj in equation (1) depend on the electron temperature Tg 
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rather than the ion temperature Ti, such difference between Ti and will affect the results 
obtained in this paper. We will address this issue in near future by considering the relaxation 
processes between ions and electrons. 
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